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0^ ' We calculate the optical conductivity of superconductors in the strong-coupling limit. In this 

0^ , anomalous limit the typical energy scale is set by the coupling energy, and other energy scales such as 

the energy of the bosons mediating the attraction are negligibly small. We find a universal frequency 
dependence of the optical absorption which is dominated by bound states and differs significantly 
^ , from the weak coupling results. A comparison with absorption spectra of superconductors with 

^0 ' enhanced electron-phonon coupfing shows that typical features of the strong-coupling limit are 

already present at intermediate coupling. 
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I. INTRODUCTION 
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According to the traditional pairing theory of superconductivity a high transition temperature may have two origins. 
The BCS-formula 

O 



ksT, w nLUoexpi-l/NpY) (1) 



suggests that the transition temperature is large if either the typical energy of the attractive pairing interaction, htoo, 
or the strength of the pairing interaction, NpV, is large. Both possibilities were taken into consideration as the origin 
of the high transition temperatures observed in cuprate superconductors. Most theories of high-Tc superconductivity 
■ invoke the first mechanism, i.e. a large typical energy of the bosons mediating the pairing interaction and a small or 
^ ' intermediate coupling constant. These systems are well described by the weak-coupling theory of BCS. Alternatively, 
O it has been suggested that the high transition temperatures and other anomalies found in the cuprate superconductjors 
J-^ ' are due to an unusually strong coupling of conduction electrons to the bosons mediating tbe attractive interactionfifl. 
^ . Such systems are the subject of Eliashberg's strong-coupling theory of superconductivitycl'El. The weak-coupling BCS 
theory is obtained from Eliashberg's theory in the limit kBTc/Tiuj^ — > 0. The opposite limit, i.e. kBTc/hujQ — + oo, is 
^ , called "strong-coupling limit" , and will be discussed in this article. 

' The strong-coupling limit of the Eliashberg theory has been of recent interestl9~El in the context of the ongoing 
high-Tc debate. Indeed, EUashberg's strong-coupling theory is among the most widely discussed candidates for a 
proper approach to high-Tc superconductivity. -This includes d-wave pairing induced by a coupling of conductiaa 
electrons to antiferromagnetic spin fluctuationso, as well as s-wave pairing due to a phonon mediated attractionll^. 
In both cases, the strong coupling of the electrons to bosonic excitations leads to a very short inelastic lifetime, 
h/r w ksTc, and one of the conditions of the weak-coupling theory {H/t <C ksTc) is violated. A substantial fraction 
of the bosons have an energy smaller than ksTc, and one expects typical features of the strong-coupling limit to 
be observable in high-Tc superconductors. With this perspective, the investigation of the formal strong-coupling 
limit, kBTc/hujQ —>■ oo, becomes more than just an interesting intellectual exercise. Indeed, we can hope to gain 
from this limit some new physical insights into superconductivity in the standard strong-coupling regime, /csTc ^ huQ. 
Superconducting properties are not universal in this regime and depend on a number of material parameters. On the 
other hand, the weak-coupling and the strong-coupling limits are universal, and a thorough understanding of both 
universal limits will be useful for understanding the anomalies observed in the intermediate cases. Moreover, although 
no strong-coupling superconductor with fc^Tc > hujQ has yet been found, one can envisage that some new physical 
systems will be discovered in the future which fall in the strong-coup Jiug .UHjiit. There are indications that important 
phonons have energies of the order of 2nkBTc in Ceo-superconductoraldO'ES, and that Cgo is on the borderline to the 
very strong coupling regime. 
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In the following we will assume, for convenience, that the attractive interaction is mediated predominantly by 
phonons. The universality of the strong-coupling limit implies that details of the phonon spectrum or the electron- 
phonon coupling do not matter, so that we can let the phonon spectrum collapse to an Einstein spectrum at an average 
phonon frequency =< >^^'^. The formal strong-coupling limit then amounts to letting go to zero. Keeping the 
electron-phonon coupling parameter. A, finite in this limit would lead to a vanishing Tc, and thus to a non-interesting 
fixed point. The strong-coupling fixed point of interest is obtained by having the coupling parameter A go to infinity 
in a way such that the energy X^^'^hfl ^ stays fixed. This limit can also be formulated in an alternative way. One 
lets the phonon spring constants go to zero (ideal softening of the phonon frequencies) , and keeps the ionic mass and 
the McMillan-Hopfield interaction parameter, rj = Np < P >jlJ'tj constant (constant electron-ion coupling). There 
is only one characteristic energy left in the strong-coupling limit, namely 



E = n^^^^^.x^^^nn, (2) 

where Np is the density of states at the Fermi energy, and M the average mass per site. The BCS fot«i**ia (^ for 
the transition temperature of weak-coupling superconductors is replaced in the strong-coupling limit byllj'tZI 



ksT, ^ 0.1827?^^ ^^ ^ . (3) 

For convenience we will take E as our energy unit. In these units the limiting critical temperature, ksTc, is equal to 
0.1827 and the energy gan .equal to 1.16Q. 

The basic new featureETEl of the strong-coupling limit is a substantially altered excitation spectrum as compared 
to the BCS weak-coupling spectrum. At zero temperature, the excitation spectrum is discrete and the density of 
states is an infinite set of Dirac delta functions. The positions and the weights of these peaks have recently been 
calculatedil. Physically they correspond to bound states. The lowest of these excitation energies determines the gap 
of the superconductor in the strong-coupling limit, and one obtains 2Aq/ ksTc = 12.7. At a finite temperature the 
excited states have a finite lifetime due to the interaction with thermally activated phonons, and the energy levels get 
broadened. The positions of the peak maxima remain essentially unchanged. 

These anomalies in the spectrum lead to various anomalies in observable properties. A number of physical properties, 
in particular thermodynamic properties including the upper oritical field and the penetration depth, have already 
been studied for superconductors in the strong-coupling limitQ. On the other hand, we know of no calculation of 
any dynamical properties in this limit. A dynamical quantity of particular interest is the infrared conductivity. The 
real part of the conductivity of strong-coupling superconductors shows, in general, a quite complicated frequency 
dependence which reflects the rich variety of different absorption channels. An infrared photon can break a Cooper 
pair or accelerate a thermally excited electron. As a consequence of the electron-phonon coimling, both processes 
may be accompanied by the absorption or emission of one or several phonons (Holstein effect )l3. It would be useful 
to distinguish the more or less universal features in the infrared spectrum (if they exist) from those which depend on 
material parameters such as the phonon spectrum. The strong-coupling limit is interesting in this respect because it 
is independent of the phonon spectrum, as was mentioned above. Therefore, the only structures which survive in this 
limit are the universal strong-coupling features. 

We calculate in this paper the infrared absorption for a superconductor in the strong-coupling limit, and discuss 
the anomalies found in this limit, as well as their relevance for superconductors with a strong electron-phonon cou- 
pling. The theoretical background of our calculations is presented in Sec. II. Our results are discussed in Sec. Ill 
which also contains the comparison of the strong-coupling limit with numerical results for the infrared absorption in 
superconductors with a strong electron-phonon coupling. Our conclusions are presented in Sec. IV. 



II. CONDUCTIVITY IN THE STRONG-COUPLING LIMIT 



We_calculate the conductivity in the strong-coupling limit by standard strong-coupling theory, as developed by 
NanxJ,. Sihs-W .andjiSwiharto, and as used, more recently, in the context of high-Tc superconductors by several 
authorsBMEl'BEl. This theory gives the conductivity in the frequency|-4p,nge uj <^ Ep/h, which|-cpvers the range 
of interest, uj « Ao/?i. We use in the following the notation of Lee et al.Eil, and Rainer and Saul£3. For isotropic 
electron-phonon scattering one can express the conductivity in terms of the ^-integrated Green's functions and the one- 
particle self-energies of the superconductor in equilibrium. The retarded (superscript R) and advanced (superscript 
A) self-energies are conveniently written in terms of the renormalized excitation energy, e^'^{£), and the renormalized 
gap function, A^^'*(e): 
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The ^-integrated Green's function is then given by 



(e) = -vr- 



(e)f3-A^-^(e)(zf2) ^ _^ 6f3 -A^-^(e)(^f2) 



where we followed a standard notation of strong-coupling theory, and introduced the gap function A(e) 



A«^^(e) = A^*'^(e)^ 



R,A, 



The conductivity a{uj) is, in the long- wavelength limit {q — 0), given by the following integral: 



cr(w) 



2Dhu 



+ OC 



del tanh 



1 



tanh 



2kBT c"{€,uj) 
e+ 1 



e_e+ + A^A^ 



(A«)2 - (6_)^/(A?)2 - (e+)2 



2kBT c^(e,w) 



e_e+ -I- A^A^ 



1 



tanh — tanh ■ 



1 



2kBT 2kBT J c'^(e,w) 

where we have used the abbreviations 



A^A 



€± — t ± ?iijj/2. 



(4) 



(5) 



(6) 



(7) 



(8) 
(9) 



and 



:-(6,c.) = .\/A«(e+)2 - g^.(e+)2 + .JA-4(e_)2 _ 



(10) 



(11) 



Np is the density of states in the normal state, vp the Fermi velocity, and D the dimension of the super-eonductor. 

The self-energies and Green's functions have been calculated in the strong-coupling limit by Combescota. We adopt 
the notation of this reference and introduce, for convenience, a complex angle ip{e) defined by 



sin (/3(e) 



A«(e) 



(12) 



We can now write the conductivity (^ in terms oiipie). We first note that c^'"^'°(e) are constants in the strong-coupling 
limit, and equal to the inverse lifetime, 



(13) 



where the last equality follows from eq. (||). The inverse lifetime is linear in T and diverges in the strong-coupling 
limit. In this way the strong-coupling limit is completely analogous to the dirty limit for a superconducting alloy. 
Hence, we follow the usual procedures of the dirty limit theory, and combine the lifetime with the prefactors of the 
integral in (^ to obtain the conductivity in the normal state. 



(TJV = 



2NFe^vlT 
D 



(14) 
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The remaining terms in the integral in (|^) can be expressed in terms of by using the relations 

A«(e) _ 1 



and 



VA«(e)2 - £2 



tan(p(e) , 



tan(/3*(e) , 



v/A^(e)2 - £2 cos (/7(e) 
A-^(£) _ 1 



cos (e) 



The relations in (16) for the advanced functions follow from (^5|) by the fundamental symmetryE2l 

A^(£) = (A^(£))*. 

By making use of (|l|), (|l|), (|l^, ([^, and the symmetries 

A^'^(e) = (A^'^(-£))* . 
one finds the following expression for the conductivity: 



^ = -1 / rf£ 

(Jn Twj 



A{e^ bj) tanh 



2kBT 



A*i 



I tanh 



2kBT 



+B{t,Lj) (tanh 



2kBT 



tanh ■ 



2kBT 



where 



and 



2A{e,uj) = 1 + tan(p(£_)tan(/3(e-)_) + 



cos (/3(£_) cos (y5(£+) ' 



2B{e,uj) = 1 + tan(p*(£_)tani^(£+)) + 



COS(^*(£_) COS(p(£+) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



Formula (|T|) gives us the complex conductivity of a superconductor in the strong-coupling limit, and will be used in 
the following section. Of special interest here is the real part of cr, which determines the contribution of the conduction 
electrons to the optical absorption. By taking the real part of eq. (M) we find: 



5Re 



de 



, £+ , £ 

tanh — - — — — tanh 



2kBT 



2kBT 



X 3mtaniy9(e_) 5r7itan(^(£+) 



cos(/3(£_) cos(/5(e+) 
Eq. ( p^ ) is the basic formula for the following discussions. 



(22) 



III. INFRARED ABSORPTION 



We first calculate the optical absorption at zero temperature from the zero temperature limit of eq. (^2|). It gives 
us the absorption in terms of the function (p(e) defined in eq. (p^. It was shown recentlyu that ip{e) is real, positive, 
and increases monotonically from zero to infinity at positive energies. This means that tan(/? and 1/ cosip are real at 
almost all energies except for a set of isolated points where tan93(£) and 1/ cosiy9(£) have poles, and consequently a 
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(5-function imaginary part. These poles determine the discrete quasiparticle energies, which are discussed in detail in 
Ref. ^. It is shown there that the quasiparticle density of states, 

iV(e) = 3mtan.^(e) , (23) 

is a sum of ^-functions. Explicitly, one finds the following quasiparticle density of states at T = 0: 

oo 

N{e) = TT ^ P„ [S{e ~ En) + S{e + E,,)] , (24) 

71=1 

where the E^s are the (positive) quasiparticle energies, and the corresponding residues P„ are given by P„ — l/ip'{En) 
. The -E„'s are obtained from the condition (p{En) = (2n — l)7r/2 for n — 1,2, ... A very good approximation at e > .4 
is 

^(e) « £!±li^ + 1,05 (25) 

TT 

The first three E^s and the corresponding residue arei: Aq = Ei = 1.16, Pi = .71, E2 = 3.04, P2 = 0.48, E3 = 
4.30, P3 = 0.34,. The gap function A(e), whose real part is an even function of e, is given at the energies En by the 
simple formulaQ A(i?„) = (— l)"i?„, which leads to the relation 

^ (-!)"• (26) 



smip{En) 



This result has important consequences for the absorption spectrum, as will be shown below. 
Eq. ( p^ ) can be written at T = in the form 

2 phuj/2 

d€N{€ + Twj/2)N[e - htj/2)x 

^ (27) 



Kecr(w) = (TAT— / deN{e + hLu/2)N{e-hoj/2)x 
1 



siniy9(e + huj /2) sanp(e — Tilj /2) 



One can now insert eqs. (24), (|2^), and one obtains 

^ea{uo) = + (-l)"+"+i]P™P„<5(£;,„ + E^ ~ huo) . (28) 

The physical interpretation of this result is quite clear. A photon of frequency m is absorbed by the superconductor, 
exciting two quasiparticles of energies Em and i?„. The strength of this absorption is determined by the joint density 
of states, 

N{e + huj/2.)N{e~huj/2), (29) 

and the coherence factor 

sin(^(e + ?ia;/2) sin(^(e — ?itj/2) ^'^^^ 

Naturally, the absorption spectrum is discrete, due to the discreteness of the quasiparticle density of states. The 
interesting new features are stringent selection rules for optical transitions. These rules originate from the coherence 
factor ( ^o| ) together with the relation (|2^), and produce the factor [1 + \^ (p8|). Thus, creation of two 

excitations with quantum numbers m and n of the same parity (n + m even) has a vanishing coherence factor and 
is optically forbidden. An interesting example of a forbidden process is the excitation threshold at energy 2Ao 
(m = n = 1). Actually, the same phenomenon, i.e. a vanishing coherence factor at threshold, arises in the Mattis - 
Bardeen theory for the optical absorption of weak-coupling superconductors in the dirty limit. Here, the vanishing 
coherence factors compensate for the diverging joint density of states and one finds a finite absorption above threshold 
which vanishes continuously if the threshold is approached from above. The excitation spectrum of the weak coupling 
superconductor is continuous. In the strong-coupling limit, on the other hand, the excitation spectrum is discrete, and 
the effect of the destructive coherence factor is more dramatic. The threshold for optical absorption is no longer at the 
'expected value' of 2Ao = 2.32. Instead, it corresponds to the creation of an excitation with energy Ei = Aq = 1.16 
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and one with energy E2 = 2.6Ao = 3.04. As a result, the threshold for absorption is at 3.6Ao, which is about a factor 
2 higher than what one might otherwise have anticipated. 

Since the absorption spectrum is a set of delta functions, it is not convenient to represent it directly on a graph. 
Rather, we display in Fig.l the integrated absorption, dLLi^ea{Lu). The positions of the steps give the positions of 
the absorption peaks, and the height of the steps give the weights of the peaks. The steps in this figure correspond 
to pair breaking into Ei + E2 ^ 3.6Ao, Ei + E^ = 5.6Ao, E2 + E3 = 6.3Ao, Ei + Eq = 6.9Ao, E2 + E5 = 7.9Ao, 
Ei+Es= 8.OA0, and E3 + E4 = 8.3Ao. 



3^ 
3 



CD 




FIG. 1. Integrated conductivity in the strong-coupling limit for T ^ 0. Frequencies are given in units of Ao/h, and the 
conductivity in units of on{T 0). The steps indicate the location of (5-peaks in the absorption, and the heights of the steps 
give the weights of the peaks. 



We now turn to the case of finite temperatures. The equation for the real part of the conductivity in the strong- 
coupling limit is given in eq. (|2^). At finite T the function (^(e) has an imaginary part, and the absorption spectrum 
becomes continuous. This function can ba calculated numerically by various techniques, for example by solving a 
second order nonlinear differential cquationcl. The results of this calculation are shown on Fig. 2 for three temperatures: 
T/Tc = 0.3, 0.6 and 0.9. The absorption for T/Tc = 0.3 is easily interpreted with the help of the T = result. The 
main absorption peak at an energy of huj ~ 3.5Ao corresponds to the T = absorption threshold at 3.6Ao. Its 
position is slightly shifted to lower frequencies, in agreement with the small shift of the peaks in the density of states. 
The next peak corresponds to the breaking of a pair into two excitations at Ei and E^. Its position at T /Tc — 0.3 is 
essentially unchanged from the T = position. In addition, we find at 6.3Ao the peak corresponding to pair breaking 
into E2+ E3. The peak corresponding to Ei + Eq is barely seen as a shoulder on the side of the preceding two peaks. 
Finally, the last peak in the figure comprises the T = peaks at 7.9Ao, 8.OA0 and 8.3Ao which are unresolved. 

In addition to these structures, which are in correspondence with the peaks at T = 0, the absorption shows for 
T/Tc = 0.3 a small threshold at about the standard absorption threshold 2Ao. This feature indicates that the 
coherence factors are no longer strictly zero, because the excitation energies acquire a finite width due to their finite 
lifetime. At T/Tc = 0.6 this 2Ao threshold has grown, but the dominant feature of the absorption spectrum is still 
the main peak at « 3.5Ao whose frequency is almost unchanged. On the other hand, the structures in the spectrum 
at higher frequencies have been essentially washed out, except for a small bump at w 6, corresponding to the fusion of 
the 5.6 and the 6.3 peaks. At T/Tc — 0.9, the gap is filled with states, and there is a substantial amount of absorption 
down to zero frequency due to the presence of thermal excitations. One clear feature survives in the spectrum at 
elevated temperatures, which is the peak at « 3.5Ao. This is the remainder of the T = threshold. 
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FIG. 2. 5Re(T(a;) in the strong-coupling limit as a function of the frequency. Results are shown for T/Tc = .3 (solid line), .6 
(short dashed line), and .9 (long dashed line). The conductivities are given in units of (Jn{T), and the frequencies in units of 
Ao/a. 



We see that the main features of the absorption at T = are quite robust, and survive thermal smearing. Thus, 
it is interesting to check if these features of the strong-coupHng Hmit are still present at finite values of the coupling 
constant, and even down to realistic couplings. In order to explore this possibility we have solved Eliashberg equations 

on the real frequency axis for a fixed shape of the spectral function, and several values of the coupling constant A. 
For convenience we have taken a narrow gaussian spectrum (Einstein-type spectrum) of the form: 



where VLq is a characteristic phonon frequency, and F = Since, as we have seen, the strong-coupling limit 

automatically implies the dirty limit, we have also used the dirty limit for our calculations at finite A. The reduced 
temperature is low in all the calculations (T/Tc = 0.22 for A = 1 and much lower for the higher A's). Figures 3 and 
4 show the results for 5Recr(a;) for A = 1,2,3,4,10 and 40, the ratio f2o/Ao being respectively 4.26, 1.83, 1.23, 0.95, 
0.46 and 0.183. We note that our program recovers the standard Mattis-Bardeen behavior for A ^ 0. 




(31) 
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3.0 




FIG. 3. ^ea{iu) in the dirty limit for Einstcin-typc spectra as a function of the frequency. Numerical results at very low 
temperatures are shown for A = 40 (T/Tc = 0.022; solid line), and A = 10 (T/Tc = 0.045; dashed line). The conductivities are 
given in units of crjv(r), and the frequencies in units of Ao/h. 




12.0 



FIG. 4. ^ea{u)) in the dirty limit for Einstein-type spectra as a function of the reduced frequency for A = 4 (solid line), 
A = 3 (long-dashed line), A = 2 (dashed line), and A = 1 (dot-dashed line). The BCS weak coupling limit is given, for 
compaxison, as the thin solid line. The arrows indicate the position of the Holstein edges for A = 2 (right arrow), A = 3 (middle 
arrow), and A = 4 (left arrow). The conductivities are given in units of a-jv(r), and the frequencies in units of Ao/fi. 



As can be seen in Fig. 3, our numerical calculations for large A (A = 40) show a very strong similarity with the 
strong-coupling limit (A = oo). The numerical results are found to be practically indistinguishable from the strong- 
coupling limit for A > 100. Characteristic features of the strong-coupling limit are still visible for smaller coupling 
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constants (see Figs. 3 and 4). However, they have smaUer amplitudes and are shifted upward in frequency. It is 
interesting to compare the location of these peaks with the Holstein thresholds located at hu = 2 Aq + fi-Oo (see arrows 
in Fig. 4). These are the thresholds for breaking a Cooper pair and, at the same time, emitting a phonon. We see 
that the positions of the absorption peaks have at intermediate coupling constants (A^ 2) no relation to the Holstein 
threshold. 



IV. CONCLUSION 

Our results for the infrared absorption in the limit ksTc/fuvo oo (strong-coupling limit) revealed some surprising 
features. The absorption threshold of traditional weak and strong-coupling superconductors is at 2Ao, which is the 
energy needed to break a pair into two electrons and to put them into the lowest excited quasiparticle state at energy 
Aq. In contrast to this, we found that the threshold in the strong-coupling limit is given by a different process. One 
electron of a broken pair is excited into the lowest energy state at Aq = S.SSfcsTc and the other one into the next 
higher state with an energy 2.6Ao. Accordingly, the threshold for infrared absorption is found at 3.6Ao instead of 
2Ao (= 12.7A;bTc). The disappearance of the absorption at 2Ao is a consequence of the coherence factors which lead 
to a vanishing transition amplitude for the 2Ao transition. Although the absorption spectrum is smeared out at finite 
temperatures, a strong absorption peak at hu w 3.6Ao remains visible in the infrared spectrum for temperatures 
almost up to the critical temperature. As in all pairing theories of superconductivity, one finds in the strong-coupling 
limit no absorption gap at a finite temperature. Absorption in the gap is brought about by thermally activated 
excitations. Absorption near zero frequency becomes important only near Tc- Otherwise, the dominant low-frequency 
feature at finite temperatures is the appearance of a small absorption at the standard threshold of 2Ao. This feature 
grows with increasing temperature and becomes comparable to the main absorption peak near Tc- 

Our results for the infrared absorption in the strong-coupling limit are quite different from standard absorption 
spectra for superconductors with a weak or intermediate coupling. In order to investigate precursor eff'ects of the 
strong-coupling limit at finite coupling, we have calculated the absorption spectra for dirty strong-coupling supercon- 
ductors with an Einstein-type spectrum and a series of different coupling strengths. We found that the characteristic 
structure at fiu ^ 3.6Ao, which corresponds to the main absorption peak in the strong-coupling limit, is clearly visible 
for A's down to 3. It shows up as a 20% bump over the normal state value. However, the position of this peak 
shifts upwards progressively when the coupling strength is lowered and, for A = 3, it is located at about 7.5 times the 
gap Ao, which exceeds the 3.6Ao found in the strong-coupling limit. 
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